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Abstract 

Formulas are developed for the eight basis matrices {T^ ,U^} of the 

finite dimensional (p, g)-irreducible representation of SU(3). A computer program is 
included that, given p and q, calculates the basis matrices. 
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1 Introduction 

When working with a group, it is sometimes convenient to have numerical matrix represen- 
tations. Explicit formulas for basis matrices of SU(3) have been available for decades. [I] 
However the earlier versions anyway are not written in the language of currently popular 
elementary particle physics texts. [2J Given the volume of literature since SU(3) became im- 
portant in particle physics it is likely that formulas are available in the popular language 
and have been coded into computer programs. However, the author is unaware of any. This 
article includes a computer program and uses the notation and quantities found in popular 
textbooks. 

We start from the properties of the states in a [p, q) irreducible representation (irrep) 
and obtain basis matrices that satisfy the well-known commutation relations of the algebra 
su(3). No attempt to show the equivalence with the earlier work of others is attempted. 

For SU(3), there are eight basis matrices. Three of the eight, popularly called and T^, 
represent an SU(2) subgroup and can be written with standard spin matrices once the T-spins 
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of the states have been determined. For each state, one can also determine by conventional 
methods the '3-component' of its U-spin, another SU(2) subgroup. The 3-components form 
a fourth, diagonal matrix, U^. Of the final four matrices, two are the transposes of the other 
two, leaving two, and V~^, to find by various means. 

The matrices f/"*" and were obtained, in part, by solving the commutation relations 
and hermiticity requirements. Less rigorously, some expressions were developed by trial and 
error using a numerical catalog of matrices for small {p, q) found by computer. 

The matrices here are written for basis states distinguished by {p, q, s, a, u^} where (p, q) 
labels the irrep, is the square of the T-spin, a is one component (the '3-component') of 
the T-spin, and is one component (another '3-component') of the U-spin. States first in 
line are those with the lowest T-spin s then by the 3-component of T-spin, largest a first. 
A set of 2s -|- 1 states with the same T-spin s but lower 3-components of U-spin precedes 
similar sets with larger U-spin 3-component u^. 

As a practical matter, those using the formulas are advised to check for the proper 
trace and hermiticity of the matrices, ensure the examples satisfy the commutation relations 
and compare the Casimir invariants with well-known formulas. Matrices obtained with the 
formulas have been checked by the author by computer for all {p, q) irreps with p > q and 
with dimension d = {p + l)(g -|- l)(p + q + 2)/2 less than 300. Matrices for {p,q) irreps 
with q > p can be found as the negative transpose of the (g, p) irrep matrices which can be 
calculated with the formulas here. 

A Mathematica [3] notebook is available at | http : / / my web . wit . edu/ shurtleff r / and is 
included in Sec. [101 



2 The Fundamental Rep 

Unitary 3x3 matrices with determinant equal to one form a group under matrix multiplica- 
tion; the product of two unitary matrices is unitary and the hermitian conjugate of a unitary 
matrix is its inverse and is unitary. 

A unitary 3x3 matrix D{6) can be written as a function of eight real parameters 6i, i & 
{1, 2, 8} with eight generators F*, 

D{e) =expieiF' . (1) 

It follows that the generators are traceless 3x3 hermitian matrices. 
A well-known set of F^s is [2] 
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^ /O 1\ ^ /O -A ^ /O 0\ 

F^ = - F^ = - F^ = - 1 
^\10 0/ ^\iO 0/ ^\0 10/ 




By inspection, these F*s are hermitian and traceless. 

Notice that {F^,F'^,F^} contain the Pauh 2x2 spin matrices of SU(2). Embedding 
standard matrices from SU(2) is a major simphfication in the construction of SU(3) matrices 
for other irreps. 



3 General Considerations 

The {p, q) irrep, labeled with nonnegative integers p and g, of SU(3) can be represented with 
a basis of eight d x d square matrices of dimension d = (p + l)(g + + g + 2)/2. [2] 
The matrices in ([2]) represent the (p, q) = (1, 0) irrep with d = 3 dimensions. 

Since the negative transpose of the eight basis matrices for the {p, q) irrep is a represen- 
tation for {q,p), we need only consider p > q and use the negative transpose for the q > p 
irreps. 

For a {p, q) irrep, the matrix D{6) in ([1]) is a transformation matrix acting on (i-component 
'vectors' each of whose components is proportional to a 'state'. The basis states are vectors 
each having just one nonzero component. [2] Two of the matrices, and U^, are here diagonal 
and so the basis states are eigenvectors of and with eigenvalues making the diagonal 
components of the matrices and U^. The eigenvalues are also called 'quantum numbers' 
and the basis states form a 'multiplet'. 

It is equivalent to work with the eight matrices {T"*", T~, T^, V^, V~ , U^, , U^}, where 

= F^± iF^ ; = ; = F^ ± iF^ ; (3) 

U^ = F^±tF' ; U' = --F' + y^F' . 

2 2 

These equations can be inverted to find the F*s from the Ts, U s and Vs. 

We require the to be (i) hermitian and (ii) traceless. Furthermore we insist that (iii) 
{F^, F^, F^, F^, F^} have real- valued components while (iv) {F^, F^, F"^} are pure imaginary. 

By (i), (iii), and (iv), it follows that (A) the {T+, T", T^, I/", f/+, f/", f/^} all have 
real-valued components. It also follows that 

(B) Tr=T^ and (C) t/," = f/^ and (D) Vr = V^^ ■ (4) 
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Therefore, once we find and we get U~ and V~ by taking the transpose. We get T~ 
another way. 

The eight matrix basis has 28 = 8 * 7/2 commutation relations, |1] 

[T^T^] = ±T^ ; [T+,T-]=2T^ ; (5) 

[T',U^] = tIu^ ; [T',V^] = ±^V^ ■ [f/3,T±] = TV ; (6) 

[U^,U^] = ±U^ ; [U\V^] = ±^V^ ; [T\U^]=0 ; (7) 

[T±, U^] = [T^, = ; [T^, ?7^] = ± ; [T^, V^] = ±U^ ; (8) 

[f/+,?7"] = 2f/^ ; [V+,V-] =2U^ + 2T^ ; [\/^, [/^] = ± ; [U^,V^]=0 . 

(9) 

Notice that the Ts (and f/s and Vs) satisfy the commutation relations of SU(2). Following 
convention, we build the three Ts from standard spin matrices. 

We write the Ts as the direct sum of irreps of SU(2) with each SU(2) irrep occupying 
a block along the diagonal with vanishing components outside those blocks. Standard spin 
matrices are available for these diagonal blocks, so the Ts are determined once the list of 
spins is known. The list of T-spins can be found by conventional methods. Conventional 
methods also give the values of for the basis states. [2] 

Therefore we consider T^, T^, and as given, even though some work is needed to find 
expressions for them. By (jlj), we have U~ and V~ as transposes of and V^. Thus the 
unknown matrices are and V~^. 



4 States and T-Spins 

We start with the states of the (p, q) irrep to find the list of T-spins. 

Since any two linear combinations of T^ and are diagonal and commute with one 
another, the states can be identified just as well by other quantum numbers. It is conventional 
to graph the states as in Fig. 1 where the coordinates are the T^ eigenvalue and the value 
of y = (4/3)f/^ + (2/3)T^. Many states have the same T^, Y coordinates and the points in 
Fig. 1 are labeled with the number of states at the point, the 'multiplicity'. 

The T^, Y plot is useful. The number of spaces in the top row in Fig. 1 is p = 5 and 
the number of spaces in the bottom row is g = 3. And, since [T^, Y] = 0, the action of the 
step-up and step down matrices T^ is to move from state to state on the same row, Y = 
constant. 



4 STATES AND T-SPINS 



5 



For example the row in Fig. 1 with Y = —4/3 has 16 states. Of those, 7 states have 
multiphcity 1 or more, 5 states have muhiphcity 2 or more, 3 states have muhiphcity 3 or 
more and 1 state has muhiphcity 4. One can infer the T-spins s for the row by setting the 
number of states equal to 2s + 1, so the set of T-spins for the states on that row are s = 
{3,2,1,0}. 

Looking at all the rows in Fig. 1, one finds the following collection of T-spins, doubled 
for convenience and arranged from small to large, for the (5, 3) irrep, 

2s = {0,1,1,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,7,7,8} . 

There are (5 + 1)(3 + 1) = 24 double T-spins 2s listed for the (5, 3) irrep. 

The T-spins are grouped into 'top cap', 'bottom cap' and 'middle'. For the (5,3) irrep, 
we have 

top cap: 2s = {0,1,1,2,2,2} ; middle : 2s = {3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5} , 

bottom cap : 2s = {6, 6, 6, 7, 7, 8} . 

Notice the symmetry; there is one smallest (2s = 0) and one largest (2s = 8), two next 
largest at 2s = 1 and two next smallest at 2s = 7, etc. 

One finds similar results for the general {p, q) irrep. There are {p + l)(g + 1) T-spins s. 
For the (p, q) irrep and doubled to 2s, we have 

2s = {0,1,1, l;g,...,p;p+ l,...,p + g- l,p + g- l,p + g} [p>q] , (10) 

where the semicolons separate the top cap, middle and bottom cap. For the (p, q) irrep, we 
have for the doubled T-spins 2s, 

top cap : 2s = {0,l,l,...,g — l,...,g— 1} , (11) 
2 total of q terms 

middle : 2s = {g, g, g + 1, g + 1 , ..., p- l,. . .,p- l ,p, ••■,p} , 
q+1 q+1 q+1 q+1 

bottom cap : 2s = {»+ 1, ...,»+ 1, + g — l,n + g — l,p + g} \p>q\ 

^ V " s/ " 

q 2 

For the (g,p) irrep with p > q, one finds the same T-spins, so the formulas (fTT]) do not 
apply to any {p, g) with g > p. In the case of a {p, g) irrep with q > p, one would need to 
interchange p and q, p ^ q, in (fTTl) to get the correct T-spin list. 
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Thus the number of T-spins Nt in the hst for the [p, q) irrep is 

1 g 
A^T = J]^ + (p-g + l)(g + l) + J]fc = + l + + = (p + l)(g + l) . 

1=1 k=l 

Since the number of states for each T-spin s is 2s + l, the number of states, i.e. the dimension 
d of the (p, q) irrep, is 

^ = 5Z^' + (^ + l)I]0' + l) + 5Z^(P + ^-^ + 2) = ^{p + l){q + l){p + q + 2) . 

2=1 j=q k=l 

This is the expected result [2] , and confirms the hst of T-spins. 

5 T-Matrices 

The basis matrices have d dimensions. The components are separated into arrays caUed 
'blocks' based on T-spin. The ij block has 2sj -|- 1 rows and 2tj + 1 columns, where Si is 
the ith T-spin in the list (11 01) and tj is the jth T-spin in the list. There are {p + + 1) 
T-spins for the {p, q) irrep, so its basis matrices have a {p+^){q + 1) by {p + l)(g + 1) block 
structure. 

Given the list of T-spins in Sec. HJ one can apply well-known standard spin matrices to 
obtain the T-matrices T"*", T~, and T^. Arrange the states by the T-spin list (JTOl) and (JTTi) 
and within each T-spin put the state with the largest spin component first. This is the order 
in the Pauli matrices, larger spin component first. 

The T-matrices are block diagonal. The matrices T"*", T~, and are given by [5] 

^(t),, = '"^(^i'^)'^^iWi ; ^(li),, = ^^M^<x,r , (12) 

where Sj is the jth T-spin in the list ( ITOll and 

r+(s,(T) = V(s-a)(s + (T + l) ; r-{s,a) = y/{s + a){s - a + 1) , (13) 

with cr, r G {s, s — 1, — (s — 1), —s} and 5o-t = 1 for cr = r and (5o-t = otherwise. Any 
component of the T-matrices outside these diagonal blocks vanishes. 
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6 

We assign addresses to components the same as for the T- matrices: U^-j-^ , with a G {si,Si — 
1, —Si} and r G {tj, — 1, —tj}, where Si and are the ith and jth spins in the T-spin 
hst (|T0|) . Since is diagonal, we have 

Ufij)^^=U^{t,(T)6^jS^r . (14) 

Thus u^{i, Si) is the lead diagonal element of the ii block of U^. 

The values of u^{i, a) are the eigenvalues of for the basis states and can be determined 
from the values of and Y obtained for a given {p, q) by conventional means. For example, 
for the (5,3) irrep in Fig. 1, for the state | T^,Y) = \ 4,2/3), the right-most point on the 
plot, we have u^{2A,A) = 3Y/A - = -3/2 . 

For any {p,q) irrep, the commutation relation [[/'^,T"*'] = — 2T+ implies a recursion for 
the components in any one diagonal block. One finds 



v!^{i,a -I) =u^{i,a) + - ; u^(i, a) = m=^(z, s^) + , 

where is the ith T-spin and there is no sum over i. Thus we need only the lead component 
of each block u^{i, Si) to get all the v?{i, a)s for the block. 

For the (p, q) = (5, 3) irrep, the lead components of the 24 diagonal blocks of are 

top cap: 2M^(ii, Si) = {-2,-4, -1,-6,-3,0} , 

middle: 2m^(2, s^) = {-8, -5, -2, 1, -7, -4, -1, 2, -6, -3, 0, 3} , 

bottom cap : 2u^(z, Sj) = {— 5, — 2, 1, — 4, — 1, — 3} , 

where the lead components are doubled for convenience and the states are reordered if needed 
so that states with the same T-spin s but lower lead value u^{i, s) are first. 

In general, for the {p, q) irrep, the lead components of the {p + l)(g + 1) diagonal blocks 
of are found to be 

top cap : 2u^{k, Sk) = —p + q — 2i + 3j , (16) 



where k = 1 + i{i + l)/2 + j, i e {0, 1, - 1} and j G {0, 1, - 1}, 

middle : 2u^{k, Sk) = —p — q + i + 3j , 
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where k = q{q + l)/2 + i{q + 1) + j + I, i E {0,1,2, ...,p - q} and j E {0, 1, 2, g}, 

bottom cap : 2u'^{k, Sk) = ~2q + i + 3j + 1 [p > q] , 

where k = {p+l){q+l)~{q-i){q-i+l)/2+j + l, i E {0, 1, ...,q-l} and j E {0, 1, ...,q-l-i}. 
The lead components are again doubled for display. 

For a {q,p) irrep with p > q, we have u^'{k, Sk) = —u^{k, Sk)— Sk- Because of the minus 
sign, there is also a reordering of the states with the same Sk- Instead of writing more 
expressions like (fT6i) for the case (p, q) with q > p, we use the negative transpose of the {p, q) 
irrep to get matrices for the (g, p) irrep with p > q. 

At this point, for the {p, q) irrep with p > q, we have a one-to one labeling of the states. 
Besides p and q, we have put states with lower T-spin first and for a given T-spin Sj states 
with larger 3-component a go first. This still leaves whole sets of 2s j + 1 states each that 
can be interchanged. The sets of states with the same T-spin Sj are ordered by U^, with a 
2si + 1 state set with lower lead u^{i, s) is first. Thus a state can be identified by the five 
quantities {p, q, s^, a, u^{i, a)} and the values of those quantities determines the order of the 
states. 

Thus, for p > q, we can determine four matrices T^, T^, and U^. It remains to find 
and V+. 

7 U+ and V+ 

To find the matrices and we first find that just one diagonal in each block can have 
nonzero components. Then we find which blocks can be nonzero. The nonzero blocks are 
off-diagonal, taking their places in the lessor diagonals above and below the main diagonal 
blocks. 

Inside a block. By the commutation relations [T'^, = {—1/2)U~^ and [T^, \^+] = 
(—1/2)1^+ together with the expressions flTT]) for T^, one finds that each or block 
can have but one nonzero diagonal, 

^(t-).. = 4.)(^)^-.-V2 and K+.)^^ = 4)(a)5.,.+i/2 , (17) 

where a E {sj, Sj — 1, — Sj}, r E {tj, tj — 1, —tj}, with Sj and tj the ith and jth T-spin 
in the list ffTOj) . Since a and r must differ by one half, if 2si is even in a nonzero block, then 
2tj is odd and vica versa. 

Recursions follow from the commutation relations [T^ , f/+] = and V"^] = 0. One 
finds that 

r"(t„a + l/2)n+ (a-l) = r-(s„a)n+ (a) , (18) 
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for and 

r+(s„ a - 1>+ )(a - 1) = r+(t„ a - 3/2>+ ^(a) , 

for f 

Therefore, by (fT7|) . the nonzero components of each block f/^"^^.^ and V|^.-| form just one 
diagonal in the block and, by fll8l) the components on that diagonal depend on the value of 
the lead component. There is just one unknown per block, one for U^^j^ and one for V^-y 

But, by the commutation relation [T"^, U^] = V^, the unknown quantity f^^-^ of the 
V^j^ block is proportional to the unknown u^^.^ for the U^^^ block, the factor depending on 

whether Si is greater than or less than tj. Thus we need to find one unknown quantity u^^.^ 
per block. 

Nonzero blocks of and . We now use the diagonal matrix . From [?7^,t/''"] = 
f/+ and [f/^ V+] = V+, one finds that the blocks U^^.^ and V^^.^ can be nonzero only for 



3 ± 1 ... ,1 

2 



u^{i, Si) - u^{j, tj) = —— with tj = Si±- [p>q] , (19) 



where u^{i, Si) is the lead diagonal component of the ii-block of the matrix. The values 
of 2u^{i, Si) are listed in ( fT6l) . 

Since Sj and tj must differ by ±1/2, the recursions f|T7j) and f|T8l) give the expressions 



a+1 



= (t, = s, -1/2); M+.^(a) = ^^-^^ {tj = Si + 1/2), 

(20) 

and 



4,)(^) = u+.^ (tj = - 1/2); v+.^{a) = ^ "foO = Si + l/2), 

(21) 

where p > q and u^^-j is a constant, one for each block. 

Roughly speaking, for {p,q) = (5,3), there is a nonzero diagonal zig-zagging above and 
another one below the main diagonal which adds up to a few less than 2 x 24 possibly 
nonzero blocks in f/+. Since the dimension d of the (5, 3) irrep is 120, the problem of finding 
the 2 X 120^ components of and is thus reduced to finding 40 or so values of u^ijy 
We turn now to that. 
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8 Finding the Block Unknowns 



One finds that the squares of the uj^ij-^s are rational. By solving for the u^---^ s, and then 



taking the square root to find the W^ij^s, one introduces a number of arbitrary ± signs. Here 
we simply take the positive root which seems to work well. 

The u'^^j'^ s form a {p + l){q + 1) square matrix. The matrix reflects the top cap, middle, 
bottom cap structure of the list of T-spins for the states of the {p, q) irrep, (ITU]) . One flnds 
that, for any {p, q) the nonzero blocks are located above and below the diagonal in each of 
the three regions. The W^ij-j^ matrix has the general structure, for any {p, q), 



/o 


utc 











\ 


Itc 


























udm 














1dm 


























ubc 


\0 











Ibc 


J 



(22) 



where we have utc - upper top cap, Itc - lower top cap, udm - upper diagonal middle, Idm 
- lower diagonal middle, ubc - upper bottom cap, Ibc - lower bottom cap. The structure is 
evident in Fig. 2. 

To flnd formulas for the W^ij^s in the top cap and on the top diagonal of the middle region, 
a catalog of low-valued {p, q) irreps was constructed by solving the commutation relations 
for each (p, q) separately. For example, m^^^)^ = p{q + 2)/2 was found by this method. 

Once the top cap and upper diagonal of the middle is obtained, the commutation relation 
, U^] = —T~ and the fact that V~ is the transpose of produce the following equation 



for the M^^.j s along the lower diagonal of the middle region. 



u 



u 



UJ-Q+O or-1) 



2 1 + 2 

2sj^(^"''+^+°™"^'^') 



1 



2s, + 1 



-u 



(23) 



where we choose '0' in the central two terms when j < no(g + 1) and ' — 1' otherwise and j G 
{no{q) + 1 or 0, ...,no(g) + {p — q + 2){q + 1) — 3}, with j starting with UQ^q) + 1 for g = 1 
and starting with no(g) for q >2. The quantity Sj is the jth T-spin and 



no(g) = 



q{q - I) 



+ 1 



(24) 



The terms in ( !23ll lie on a '+ '-shape in the matrix because they are in the jth row or the j'th 
column. Since expressions for the top cap and upper diagonal of the middle were previously 
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found using the catalog, expressions for the lower diagonal of the middle portion follow from 

m- 

With the top cap and middle done, the bottom cap was found by a combination of solving 
f l23|) and by deducing expressions from the catalog of low (p, q). 

The u^ij^ formulas. 

We separate the formulas for g = because with g = there is no top cap, no upper 
diagonal for the middle and no bottom cap. There is just the lower diagonal for the middle 
(1dm) which is given by 

1dm: M+.+^^^.)'=p-j + l [g = 0] , (25) 

where j G {1, The result follows from (l23l) because no(0) = 1 and the last two terms 

on the right side vanish. 

For g > 1, one finds formulas for the top cap. 



utc: 4,+,,) =^^^^^ (p+no(gi+l)-2+l)(g-no(gi+l)+«+l) [g > 1] ,(26) 



2 UQ^ qi + 

qi + 1 

where qi G {1,2, ...,g - 1}, i G {'n.o(9i), •••,'n-o(Q'i + !)}• We have 

Itc: = ^ qi{qi + '^) ^ ^ ^^^^^^ " ^^^^ " ^^^^^^ + + ' 

where qi G {1, 2, g — 1}, j G {nQ{qi), 'rio(gi + 1) — !}• In addition, 

2 1 

Itc : M+ ^) = -p{q + 2) . 
The middle region u^ij^^ are found to be 

udm: ^ g + 1 + /(z^ g) + ^ + ^ ~ ^))(^ + ^)) ' (^7) 

where z G {^0(0'), •••5 ^ — lio. + l)/2} and 

■m{i,q) = i — no{q) mod g + 1 ; /(i, g) = floor(- — 

We have 

11 + 2 I1+ ^_L^+ ^ + 2 

mm: - -i + + ^M(j_5+l+0or-l,i) ~ 2s + l"(j'J+9) 



9 RECONSTRUCTION 



12 



where we choose '0' in the central two terms when j < nQ{q + 1) and '—1' otherwise and we 
have j G {^o(q') + 1 or 0, riQ^q) + (p — g + 2)(g + 1) — 3} with j starting with UQ^q) + 1 for 
q = 1 and starting with nQ{q) for q >2. 
For the bottom cap, we have 

ubc: u+.^^^)' = (r2o(gi + l)+2-Ar + gi-l)(g-gi-no(gi + l)+Ar + 2-z) X (28) 

^ p + q - qi- no{qi + 1) + + 3 - i 
p + q - qi + 2 

where gi G {1, 2, g - 1}, i e {N - no(gi + 1) - gi + 1, - no(gi) - gi + 1}. 

: wj+^i+ij)^ = (N - qi - no{qi) + 1 - j){p - N + qi + no{qi) + j)x 

^^ p + q + j - N + qi + no(gi) + 1 
(p + g - gi + l)(p + g - gi + 2) 

where gi G {1, 2, g — 1}, j G {A^ — rio(gi + 1) + 1 — gi, — riQ^qi) — gi}. These formulas 
complete the solution. 

9 Reconstruction 

{p, g) irrep with p > q. 

The formulas are written for a (p, g) irrep with p > q. The w^^.^^s in fl25|) - fl28l) go into 
( l20l) and ( l2Ti) to make f/"*" and V"''. Then f/~ and V~ are the transpose matrices of and 
V'*'. The T^, T^, and t/^ are found in ( |T2l) . ( fT4l) and ( |T6l) . One now has all eight matrices 
{T±, T^, f/±, f/3, The eight F's can be found by (Q. 

{p, q) irrep with q > p. 

The formulas in the text give incorrect results when (mis) applied for q > p. For example, 
the T-spins of the {p, g) and {q,p) irreps are the same. One can see this from a T^, Y plot of 
the states hke Fig. 1. For the (p, g) irrep, the top row has p spaces and there are g spaces on 
the bottom row. The (g, p) irrep has the same rows as but in different order with g spaces on 
the top and p spaces along the bottom row. Since the T-spin list depends on the rows and 
not their order, the T-spins are the same for (p, g) and (g,p). But (ITUI) would give different 
T-spins for the two irreps. Thus we either write the correct formulas for g > p or we try 
some other way. We go another way. 
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For q > p we first calculate matrices as above for {p' , q') = {q,p), since p' = q > p = q' . 
Then the negative transpose changes the [p' , q') matrices to an irrep for [p, q). We have 

{^jfcj ^/fcj ^fcj ^jfe' ^ifc} = fcj5 kj^ ~U kj^ ~U kj] (29) 

and 

Fjk = ~F'kj = jk 5 (30) 

where (|3]) gives the and since the F'* are hermitian, the transpose of F'^ is equivalent to 
the complex conjugate of F'^. 

When identifying the states, it is easier to work with the complex conjugate, Fj^ = 

—F'^*i^, because the order of the states is not changed. The states for {p', q') are identified 

by the labels | p' ,q' , s' ,a' ,u'^) . By definition, {p',q') = {q,p)- Since {p,q) and {q,p) have the 
same list of T-spins, we have s' = s. The minus signs in the real, diagonal matrices F^^^ = 

—F'^jf^ and Fji^ = —F'^-^ give minus signs for cr and m^. We have a' = —a and u'^ = —v?. 
Relabeling the state | p\ g', s', a', u'^), we get | q,p, s, — cr, —u^). 

We must also consider the order of the states to know what state j is connected to which 
state k by the matrix element F*^. By construction, the ordering in the matrices is the 
ordering of the primed system for {p',q'). Since the T-spins are the same, the blocks have 
the same order of T-spins Sk- But the minus sign a' = —a reverses the order for the states 
with the kth T-spin Sk, so that these 2sk + 1 states have 3-components a G {—Sk, Sk}, 
i.e. the lowest 3-component first which is just the reverse of the p > q case. States with the 
same T-spins but higher values of occur first because of the minus sign in u'^ = —u^. 
Further complicating matters, is the fact that the lead T-spin 3-component changes sign s 
—s. Then the lead component of the kk block of has u^'{k,Sk) = —v?{k,Sk)— Sk- Thus 
the labeling of the states and the order of states can be determined for the {p, q) irrep with 
q> p. 
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10 Mathematica Notebook 

This Mathematica ^ notebook calculates the eight basis matrices {Tp, Tm, T3, Up, Um, U3, 
Vp, Vm) (p is + and m is — ) for the (p,q) irrep of the su(3) algebra. The choice of p and q 
is made by the user. The notebook is available at I http : / / my web . wit . edu/ shurtleff r / 1 . 

The notebook checks that the matrices {Tp,Tm,T3,Up,Um,U3,Vp,Vm} obey the 
su(3) algebra and that they give the correct Casimer operator CI. The eight basis ma- 
trices are also calculated, but are not explicitly checked in the program. 

The LaTeX was generated by Mathematica by saving the notebook using the 'Save as 
with the 'LaTeX Document' option. 

Program Description: 

Title: SUSMatrices.nb Operating systems: Any system that supports Mathematica 

7.0. Tested with Windows XP. File Size: 92 kB 

Memory Used In Session: MUIS depends on {p,q). MUIS = 7.5 Mbyte for (1, 1); MUIS = 
8.4 Mbyte for (3,2); MUIS = 14.2 Mbyte for (5,3). 

Explanation of Quantity Names: 

Tspins - list of spins, from small to large. Tspins determines the block structure of the 
matrices. The T-spin list has structure: {top cap, middle, bottom cap}. 
Y3s - the list of the first Y3s in each block. 

U3firsts - U3firsts is the list of the first value U3 in each block. U3 is a diagonal matrix. 
The states have U3 = 3Y^/A — T^/2 . When two T3 blocks have the same spin, the block 
with the smaller U3first precedes the block with the larger U3first. 

numTspins - the number of spins in the list Tspins. The components of a basis matrix are 
organized into blocks with numTspins as the number of blocks on a side. 
block - The ij block is a matrix with 2sj+l rows and 2tj+l columns where Si is the ith T-spin 
and tj is the jth T-spin 

dimREP - the number of states in the (p,q) irrep (multiplet). 

nO - number related to size of the caps of the T-spin list. Used as initial values in Sums and 
Tables. 

mod and floor - these help organize the 'middle' region of the T-spin list which has the same 
spin repeated q+1 times. Used in formulas for upc2[s,t]. 

rp and rm - standard SU(2) functions; p stands for + (plus) and m stands for — (minus). 

{Tplus, Tminus, Tthree, Uplus, Uminus, Uthree, Vplus, Vminus} are block by block for- 
mulas. 

{Tp, Tm, T3, Up, Um, U3, Vp, Vm} are the matrices sought. 
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Upa and Vpa are the Up and Vp matrices but in terms of unknowns upc[i,j] . 

upc[i,j] - unknowns. One upc[i,j] for each block. The unknowns upc[i,j] are most directly 

related to Up. 

uppertopcap and lowertopcap arc lists of formulas for the top cap portion of the upc2[i,j] 
matrix, where upc2[i,j] — upc^L^j]'^ 

upperdiagonal and lowerdiagonal are the lists of upc2[i,j] formulas for the middle portion of 

the upc2[i,j] matrix 

upperhottomcap and lowerbottomcap are the lists of upc2[i,j] formulas for the bottom cap 
portion of the upc2[i,j] matrix 

miscSUBS are miscellaneous substitutions for upc2[i,j] that are intoduced in the above lists 
and vanish, such as upc2[2,0] (there is no column). Also included are upc2[i,j] formulas 
that would give indeterminant results using the formulas in the above lists, for example, 
upc2[3,l] involves spin which causes trouble in the formulas as written. 
upc2SUBS collects the upc2 formulas for substitution. 

upcSUBS - Since the matrices depend on upc[s,t], not the squares upc2[s,t], we need a list 
of substitutions for the upc[s,t] based on the squares upc2[s,t]. 

Begin Program: 

Note: Input is in bold print, output is in hght print. Remarks: (* Remark *). 

(*Choose (p,q):*) 

p = 3-q = 2; 

If[p > q, pqOK = Tbie, pqOK = False]; 

If[pqOK, "sweet", pO = p]; If[pqOK, "sweet", qO = q]; 

If[pqOK, "sweet", p = qoj; If[pqOK, "sweet", g = poj; 

If[pqOK, Prmt["The formulas are set up for p>q - no problem"]. 
Print ["The formulas are set up for p>q, so the notebook calculates with {p,q} = ", 
{p, q}, " and then takes the negative transpose to get the matrices for {pO,qO} = ", 
{pO,qO}]] 

The formulas are set up for p>q - no problem 

nimiTspins = (p + + 1); 

dimREP = i(p + l)iq + l)(p + q + 2); 

(*Useful functions*) 

nO[q_]:=nO[g] = + 1 

mod[n_, q_]:=mod[n, q] = Mod[n — nO[g], g + 1] 

floor[n_, q_]:=floor[n, q] = Floor[(n - nO[g])/(g + 1)] 

Unprotect[Up]; 

(^Standard Spin functions : *) 

rp[8_,o-_]:=rp[s,o-] = y/(s - a)(s + a + 1); rm[s_,cr_]:=rm[s,o-] = y/(s + a)(s - a + 1) 
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(*Labels for the states (basis vectors)*) 

Tspins = 

iFlatten[{Table[82 - 1, {82, q}, {m, 82}], Table [s2, {s2, {m, q + 1}], 

Table[s2, {s2,p + l,p + g}, {to, g + 1 - 82 + p}]}]; 
Y3s = Flatten[{Table[-2(p - q) - 3(« - 1) + 6(j - 1), {i, q}, {j, i}], 

Table[Table[-2p -q + 3i + 6j, {j, 0, g}], {i, 0,p - g}], 

Table[3 + p - 4g + 3(i - 1) + 6(i - 1), (i, g}, {i, g - i + 1}]}]; 
U3fir8t8 = 

iFlatten[{Table[-(p - g) - 2{i - 1) + 3(j - 1), {i, g}, {j, i}], 

Table[Table[(-p -q + i + 3j), {j, 0, g}], {i, 0,p - g}], 

Table[(l - 2g + (i - 1) + 3(j - 1)), {i, g}, {j, g - i + 1}]}]; 
(^Standard spin matrices : *) 
Tplu8[8_]:= 

Tplu8[s] = Table[If[a == al + 1, rp[s, al], 0, ah], {a, s, -s, -1}, {al, s, -s, -1}] 
Tminus[8-]:— 

Tminusfs] = Table[If [cr == al — 1, nn[s, crl], 0, ah], {a, s, — s, —1}, {<t1, s, — s, —1}] 

Tthree[s_]:=Tthree[s] = Table[If[(r == al, a, 0, ah], {a, s, — s, —1}, {al, s, — s, —1}] 
{Tplus[l]//MatrixFonn, Tininus[l]//MatrixForm, Tthree[l]//MatrixForm} 

/0\/2 0\/0 0\/lOO\] 
V^ , \/2 , 00 \ > 

yoo o/voa/2o/voo-i/^ 

(* Tp, Tm, T3 : *) 

Tp = Table[0, {i, dimREP}, {j, dimREP}]; 
For[s = 1, s < numTspins, 

Table[Tp[[i + Suin[2Tspiiis[[si]] + 1, {si, s-l}],j + Sum[2Tspins[[si]] + 1, {si, s - 1}]]] - 
Tplus[Tspiiis[[s]]][[i,i]], {i, 2Tspiiis[[s]] + 1}, {j, 2Tspins[[s]] + 1}]]; 
Tm = Table[0, {t, dimREP}, {j, dimREP}]; 
For[s = 1, s < numTspins, 

Table[Tm[[i + Sum[2T8pin8[[8i]] + 1, {8i, s - l}],i + Sum[2Tspin8[[8i]] + 1, {si, s - 1}]]] = 
Tminus[Tspin8[[s]]][[i, j]], {i, 2Tspin8[[s]] + 1}, {j, 2Tspiiis[[s]] + 1}]]; 
T3 = Table[0, (i, dimREP}, {j, dimREP}]; 
For[s = 1, s < numTspins, 

Table[T3[[t + Sum[2Tspin8[[si]] + 1, {si, s - + Sum[2T8pin8[[si]] + 1, {si, s - 1}]]] = 
Tthree[Tspins[[s]]][[i, j]],{i,2Tspins[[s]] + 1}, {j, 2Tspins[[s]] + 1}]]; 
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(* U3: *) 

U3 = Table[0, {i, dimREP}, {j, dimREP}]; 
For[s = 1, a < numTspins, 

Table[U3[[« + Sum[2Tspms[[si]] + 1, {si, s - 1}], i + Siim[2Tspins[[si]] + 1, {si, s - 1}]]] = 
U3fir8ts[[s]] + i(i - 1), {i, 2T8pins[[s]] + 1}]] 
(*Block unknowns *) 
uppertopcap — 
Table[ 

Table [{upc2[i, t + ql] ^ (no[qi+i]-i) + (i,o[ql + 1] - i) + l)(g - (nO[ql + 1] - i) + 1) } , 

{i,nO[ql],nO[ql + l]}],{ql,g-l}]; 
lowertopcap = 

Table [Table [{upc2[j + ql + IJ] ^ ^-°%+^ (p - {j - nO[ql]))(g + {j - nO[ql]) + 2)} , 

{j,nO[ql],nO[ql + 1] - 1}], {ql, g - 1}]; 
nOrowlowertopcap = Table['Ikble[upc2[nO[ql], m] —* 0, {m, nO[ql]}], {ql, g + 1}]; 
upperdiagonal = 

Table [{upc2[i, i + q\^ (-g^^(p + q + l- mod[i, q])(l + mod[i, q])) } , 

|i,nO[g], numTspins — ; 

lowerdiagonal[p_, q-]:= 
Flatten[ {Table[ 

{upc2b- + q + lj]^-l + upc2|j,i -q + < nO[q + 1], 0, -1, jali]]+ 
2^^^upc2[j -q+l + li\j< nO[g + 1], 0, -1, jah],^]- 

2Tspii[b]]+l ^P^^t7,i + g]}> 

{j, nO[q] + If[q == 1, 1, 0, ahk], nO[g] + (p - q + 2){q + 1) - 3}]}]/;q > 1 
lowerdiagonal[p_, 0]:=Flatten[{Table[{upc2[; + lj]^p-j + 1}, {j, l,p}]}] 
upperbottomcap = 

Table[ Table[ {upc2[i, i + ql] ^ 

((nO[ql + 1] + i - numTspins + ql - l){q - ql - nO[ql + 1] + numTspins + 2-i) 
{p + q-ql- nO[ql + 1] + numTspins + 3 - i))/{p + g - qH- 2)}, 
{i, numTspins — nO[ql + 1] — ql + 1, niunTspins — nO[ql] — ql + 1}], {ql, g — 1}]; 
lowerbottomcap 

Table [ Table[ {upc2|> + ql + 1, j] -»• 

((+numTspins — ql — nO[ql] + 1 — i)(p — numTspins + ql + nO[ql] + j) 

(j> + q + j- numTspins + ql + nO[ql] + l))/((p + g - ql + l)(p + g - ql + 2))}, 
{j, numTspins — nO[ql + 1] + 1 — ql, numTspins — nO[ql] — ql}], {ql, g — 1}]; 
miscSUBS = {upc2[a_, 0] ^ 0, If [g 1, upc2[3, 1] i(p)(g + 2), upc2[0, b_] ^ O] } ; 
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upc2SUBS = 

Flatten[{uppertopcap, lowertopcap, upperdiagonal, lowerdiagonal[p, q], 
nOrowlowertopcap, upparbottomcap, lowerbottomcap, miscSUBS}]; 
(* Up:*) 

uplus[8_,t_,<7_]:= 

uplusfs, t, cr] = If > t, up[s, t]\/s — a, If < t, up[s, t] \J ^t^a ^ ah\j^ , huhuj 
Uplus[s_,t_]:= 

Upliis[s, t] = Tablepf[o- == p — 1/2, upliis[s, t, cr], 0, ah], {a, s, — s, —1}, {p, i, — i, — 1}] 

Upa = Table[0, {«, dimREP}, {j, dimREP}]; 
For[s = 1, s < numTspins, s++, 
For[t = l,t < numTspins, t++, 

If[(Tspins[[t]] == Tspins[[s]] + l/2)&&(U3firsts[[t]] - U3firsts[[s]] == -1), 
Table[Upa[[i + Smn[2Tspins[[si]] + 1, {si, s - 1}], 
3 + Sum[2Tspins[[ti]] + 1, {ti, t - 1}]]] - 
Uplu8[Tspin8[[s]],Tspins[[t]]][[i,i]]/. 

{up[Tspins[[s]],T8pins[[t]]] ^ upc[s,*]},{i,2Tspins[[s]] + 1}, 
0,2Tspins[[t]] + l}]]]]; 
For[s = 1, s < numTspins, 
For[t = l,t < numTspins, t++, 

If[(Tspins[[t]] == Tspins[[s]] - l/2)&&(U3firsts[[s]] - U3firsts[[t]] 1/2), 
Table[Upa[[i + Sum[2Tspins[[si]] + 1, {si, s - 1}], 
3 + Sum[2Tspins[[ti]] + 1, {ti, t - 1}]]] = 
Uplus[Tspins[[s]] , Tspins[[t]]] [[i, j]]/. 

{up[Tspins[[s]],Tspins[[t]]] upc[s,t]},{i,2Tspins[[s]] + 1}, 
{i,2T8pins[[t]] + l}]]]]; 
VarUpa = Variables[Upa]; 
upcSUBS - 
Flatten[ {Table[ VarUpa[[i]] 

((VarUpa[[t]]^/.Flatten [Table [upc[i, upc2[i, j], {i, numTspins}, 
{i, numT8pin8}]])//.upc2SUBS), {i, Length[VarUpa]}]}]; 
Up = Upa//.upcSUBS; 
(*Vp*) 

vplusfs-, t_,<T_]:= 



vplus[s, t, cr] = If s > t, up[s, t\-\/s + cr. If s < t, — up[s, t\ J ^^r^, ahv , huhv 



Vplus[s_,t_]:= 

Vplus[s, t] — Table[If[cr == p + 1/2, vplus[s, t, a], 0, ah], {cr, s, —s, —1}, {p, — t, — 1}] 
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Vpa = Table[0, {i, dimREP}, {j, dimREP}]; 
For[s = 1, s < numTspins, 
For[t = 1, t < numTspins, t++, 

If[(Tspins[[t]] == Tspins[[s]] + l/2)&&(U3firsts[[t]] - U3firsts[[s]] == -1), 
Table[Vpa[[i + Sum[2T8pins[[si]] + 1, {si, s-1}], 
j + Sum[2Tspins[[ti]] + 1, {ti, t - 1}]]] = 
Vplus[Tspins[[s]],Tspins[[t]]][[i,j]]/. 

{up[Tspins[[s]],Tspins[[t]]] ^ upc[s,t]},{i,2Tspins[[s]] + 1}, 
{j,2Tspins[M] + l}]]]]; 
For[s = 1, s < numTspins, s++, 
For[t =l,t < nvunTspins, 

If[(Tspins[[t]] == Tspins[[s]] - l/2)&&(U3firsts[[s]] - U3firsts[[t]] 1/2), 
Table[Vpa[[i + Sum[2Tspins[[si]] + 1, {si, s - 1}], 
j + Sum[2Tspins[[ti]] + 1, {ti,t - 1}]]] = 
Vplus[Tspins[[s]],Tspins[[t]]][[i,i]]/. 

{up[Tspins[[4], Tspins[[t]]] upcis,t]}, {i, 2T8pins[[s]] + 1}, 

{i,2Tspin8[[t]] + l}]]]]; 
Vp = Vpa//.upcSUBS; 
(*Um and Vm *) 

Urn = Table[0, {i, dimREP}, {j, dimREP}]; 

Table[Um[[i,i]] - Up[b',i]],{i, dimREP}, {i, dimREP}]; 
Vm = Table[0, {i, dimREP}, {j, dimREP}]; 

Table[Vm[[t, j]] = Vp[L7,i]], {i, dimREP}, {j, dimREP}]; 
(*For the case q> p:*) 

If[pqOK, "sweet", Tp = -Transpose [Tp]]; If[pqOK, "sweet", Tm = -UransposefTm]]; 
If[pqOK, "sweet", T3 = -TraDspo8e[T3]]; If[pqOK, "sweet". Up = -TVanspose[Up]]; 
If[pqOK, "sweet", Um = -Transpose [Um]]; If[pqOK, "sweet", U3 = -'IVanspose[U3]]; 
If[pqOK, "sweet", Vp = -IVansposefVp]]; If[pqOK, "sweet", Vm = -Ttanspose[Vm]]; 
If[pqOK, "sweet ",p = pO]; If[pqOK, "sweet", g = qO]; 
(♦Check:*) 

Print["Check, for {p,q} = ",{p,g}, " . dimREP = ", dimREP, " {0} is good: ", 

Union[ Flatten [SimpUfy [{T3.Tp - Tp.T3 - Tp, T3.Tm - Tm.T3 - (-Tm), 
T3.Up - Up.T3 - (-|Up) , T3.Um - Um.T3 - (+iUm) , T3.U3 - U3.T3, 
T3.Vp - Vp.T3 - iVp, T3.Vm - Vm.T3 - (-|Vm) , Tp.Tm - Tm.Tp - 2T3, 
Tp.Up - Up.Tp - Vp, Tp.Um - Um.Tp, Tp.U3 - U3.Tp - (+|Tp) , Tp.Vp - Vp.Tp, 
Tp.Vm - Vm.Tp - (-Um), Tm.Up - Up.Tm, Tm.Um - Um.Tm - (-Vm), 
Tm.U3 - U3.Tm - (-|Tm) , Tm.Vp - Vp.Tm - Up, Tm.Vm - Vm.Tm, 
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U3.Up - Up.U3 - (+Up), U3.Um - Um.US - (-Urn), U3.Vp - Vp.US - (+|Vp) , 
U3.Vm - Vm.U3 - (-iVm) , Up.Um - Um.Up - (2U3), Up.Vp - Vp.Up, 
Up.Vm - Vm.Up - (+Tm), Um.Vp - Vp.Um - (-Tp), Um.Vm - Vm.Um, 
Vp.Vm - Vm.Vp - (2U3 + 2T3), 

|(Tp.Tm + Tm.Tp) + T3.T3 + i(Vp.Vm + Vm.Vp) + i(Up.Um + Um.Up)+ 
i(2U3 + T3).(2U3 + T3) - ((^^^2±s£ + p + g) IdentityMatrix[diinREP]) }]]]] 

Check, for {p,q} = {3, 2} . dimREP = 42 {0} is good: {0} 

Prmt["The notebook found basis matrices {Tp,Tm,T3,Up,Um,U3,Vp,Vm} for the ", 

{p, q}, "-irrep."] 

The notebook found basis matrices {Tp,Tm,T3,Up,Um,U3,Vp,Vm} for the {3, 2}-irrep. 
(*The F* basis matrices*) 

Fi = {l(Tp + Tm), f (Tp - Tm), T3, i(Vp + Vm), f (Vp - Vm), i(Up + Um), 
f (Up - Um), y|U3 + ^T3} ; 

MaxMemory Used [] 

8420680 
End Program. 
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A Exercises and Problems 

This Appendix would likely be jettisoned before the article would be published. Some 
answers may be found by downloading and viewing the LaTeX source of this Appendix. 

1. Rewrite the commutation relations ([6]) - ([9]) with T'^ = 2T^ and If^ = 2U^. See M. F. 
O'Reilly, Ref. 4. 

2. A numerical catalog of matrices for various {p, q) irreps gives the following table of values 
for the square of one of the f/^ block unknowns W^ij) ■ Find a formula in terms of p and q 
for the unknown ut-^. Which unknown is it, i.e. what are i and j? 





(1,1) 


(2,1) 


(3,1) 


(4,1) 


(5,1) 




+ 2 


3/2 


3 


9/2 


6 


15/2 
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(5,2) 


(6,2) 
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(3,3) 


(4,3) 


(5,3) 
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%.^) 






15/2 


10 


25/2 












(4,4) 


(5,4) 












12 


15 





3. (a) From Sec. [8|, determine the u^..-^ for the (2, 0) irrep. Then find the eight matrices 
{T±, T^, f/±, f/3, V^}. (b) Repeat for the (1, 1) irrep. 

4. Refer to Fig. 1. Find the {T^,Y} coordinates (a) of the state with T-spin s = 0; (b) 
of the two sets of two states each with T-spin s = 1/2; (c) of the nine states with T-spin 
s = 4; (d) of the two sets of eight states each with T-spin s = 7/2. 

5. To obtain the m^^.^ for the (5, 3) irrep, the positive square roots of u^ij^ were chosen. More 
generally, some or all of these square roots could have been negative. Consider all sets of 
positive and negative roots. Would all choices of the it signs give matrices that satisfy the 
commutation relations? 

6. Solve the commutation relations for the top cap and upper diagonal of the middle. Do 
the results agree with those displayed in Sec. [HP 
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Figure 1: Multiplicity of States for the (5,3) Irrep. The number of states with quantum 
numbers {T^,Y} where is the '3-component' of T-spin and the 'hypercharge' Y is given 
by Y = (4/3)[/3+ (2/3)T3, with is the '3-component' of U-spin. There are d = 120 states 
total. 
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Figure 2: The Matrix of Block Unknowns u^ij-^ for the (5, 3) Irrep. The components of 
the d X d basis matrices are sorted into 'blocks' by T-spin and U-spin, where T- and U-spin 
are SU(2) subgroups of SU(3). For the (5,3) irrep, there are 24 T-spins, including repeats, 
making a 24 x 24 array of blocks. Recursions and other constraints reduce the number of 
unknowns from d'^ — 120^ to about 40 for the (5. 3) irrep shown here. Note that the nonzero 
W^ij^^ form upper and lower diagonals in the middle region, enclosed by top and bottom 
'caps' with staggered upper and lower diagonals of their own. 



